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ABSTRACT 

We use an automated Voigt-profile fitting procedure to extract statistical properties 



> 

in 

I of the Lya forest in a numerical simulation of an = 1, cold dark matter (CDM) 

I universe. Our analysis method is similar to that used in most observational studies of 

' the forest, and we compare the simulations to recently published results derived from 

CL|' Keck HIRES spectra. With the Voigt-profile decomposition analysis, the simulation 

Q ■ reproduces the large number of weak lines (A'^hi ^ lO^^cm"^) found in the HIRES 

spectra. The column density distribution evolves significantly between z = 3 and 
, z = 2, with the number of lines at fixed column density dropping by a factor ~ 1.6 in 

the range where line blending is not severe. At z = 3, the 6-parameter distribution 
^ \ has a median of 35 kms~^ and a dispersion of 20 kms~^, in reasonable agreement 

5h ' with the observed values. The comparison between our new analysis and recent data 

strengthens earlier claims that the Lya forest arises naturally in hierarchical structure 
formation as photoionized gas falls into dark matter potential wells. However, there 
are two statistically signficant discrepancies between the simulated forest and the 
HIRES results: the model produces too many lines at z = 3 by a factor ~ 1.5 — 2, and 
it produces more narrow lines (6 < 20 kms~^) than are seen in the data. The first 
result is sensitive to our adopted normalization of the mean Lya optical depth, and 
the second is sensitive to our assumption that helium reionization has not significantly 
raised gas temperatures at z = 3. It is therefore too early to say whether these 
discrepancies indicate a fundamental problem with the high-redshift structure of the 
$7 = 1 CDM model or reflect errors of detail in our modeling of the gas distribution or 
the observational procedures. 
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1. Introduction 



Absorption lines in quasar spectra, especially the "forest" of Lya lines produced by 
concentrations of neutral hydrogen, are uniquely suited for probing structure formation in the 
high-redshift universe. The absorbers trace relatively pristine baryonic material over a wide 
range of redshifts, densities, temperatures, and ionization states. Recent advances in computer 
technology have expanded our ability to predict the conditions of the absorbing gas, while 



high-precision observations made using the HIRES spectrograph (Vogt, et al. 1994) on the 10m 
Keck telescope have quantified the statistics of the low column density absorbers to unprecedented 
accuracy (e.g., Hu et al. 1995, hereafter HKCSR| ). These data provide stringent constraints on 
theories of structure formation and the state of the intergalactic medium (IGM) at high redshifts. 
However, a detailed confrontation between theory and observations requires that simulations 
and observed quasar spectra be analyzed using similar techniques. In this paper we apply an 
automated Voigt-profile fitting algorithm to Lya spectra from a simulation of the cold dark matter 
(CDM) scenario ([Peebles 1982| ; [Blumenthal et al. 1984|) . We compare the statistics of the resulting 



line population to those derived by HKCSR from HIRES spectra. 



Recent cosmological simulations that incorporate gas dynamics, radiative cooling, and 
photoionization reproduce many of the observed features of quasar absorption spectra, suggesting 
that the Lya forest arises as a natural consequence of hierarchical structure formation in a 
universe with a photoionizing UV background ( |Cen et al. 1994 ; Zhang, Anninos, & Norman 1995 



Hernquist et al. 1996, hereafter |HKWM| ; [Katz et al. 1996b| ; [Miralda-Escude et al. 199^ ). Most 
of the low column density lines are produced by structures of moderate overdensity that are far 
from dynamical or thermal equilibrium, blurring the traditional distinction between the Lya forest 



and Gunn-Peterson (1965) absorption from a smooth IGM. HKWM used a simple flux threshold 
algorithm to identify lines in their simulated spectra, defining any region with transmitted flux 
continuously below a specified threshold as a single line. They showed that their simulation 
of a CDM universe reproduced the observed abundance of absorption systems as determined 
by Petitjean et al. (1993, hereafter |PWRCL| ) quite weh over most of the column density range 



lO^'^cm-^ 



< iVni < lO^'^cm with a significant discrepancy for A'hi ~ 10 



The traditional technique for identifying and characterizing quasar absorption lines is to fit 
spectra by a superposition of Voigt profiles. The HIRES spectra have very high signal-to-noise 
ratio and resolution, and most of the lines found in this way are weak absorbers with column 
densities A'hi < lO^^cm"^. The threshold and Voigt-profile procedures behave very differently in 
this regime, since a single feature identified by the threshold method will often be decomposed into 
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a blend of weaker lines when it is modeled as a superposition of Voigt profiles. In order to compare 
to published line population statistics from HIRES data, therefore, it is essential to analyze the 
simulated spectra by Voigt-profile decomposition. 

The physical model implicit in the decomposition technique is that of a collection of discrete, 
compact clouds, each characterized by a single temperature (or at least by a single velocity 
dispersion, which could include contributions from thermal motion and from Gaussian-distributed 
"turbulent" velocities). The simulations undermine this physical picture because the absorbing 
systems merge continuously into a smoothly fluctuating background, often contain gas at a 
range of temperatures, and are usually broadened in frequency space by coherent velocity 
flows that do not resemble Gaussian turbulence. Nonetheless, any spectrum can be described 
phenomenologically by a superposition of Voigt-profile lines, with the number of components 
increasing as the signal-to- noise ratio improves and more subtle features must be matched. The 
distributions of fitted column densities and 6-parameters provide a useful statistical basis for 
comparing simulations and observations, and this is the approach that we adopt in this paper. We 
will discuss the correspondence between the parameters of the Voigt-profile components and the 
physical state of the absorbing gas elsewhere (Dave et al., in preparation). 



2. Simulation and Artificial Spectra 

The simulation analyzed here is the same as that of HKWM: a CDM universe with fl = 1, 
Hq = 50 km s~^Mpc~^, baryon fraction Q}, = 0.05, and a periodic simulation cube 22.222 comoving 
Mpc on a side containing 64^ gas particles and 64^ dark matter particles, with individual particle 
masses of 1.45 x 10^ and 2.8 x 10^ Mq, respectively. The power spectrum is normalized to 
CTg = 0.7, roughly the value required to reproduce observed galaxy cluster masses ([Bahcall & 



Cen 1992 ; White, Efstathiou, fc Frenk 1993 ); we will consider COBE-normalized CDM models 
with = 1 and < 1 in future work. We use the N-body -|- smoothed-particle hydrodynamics 
code TreeSPH ( [Hernquist fc Katz 198S| ) adapted for cosmological simulations (Katz, Weinberg, &: 



Hernquist 1996, hereafter KWH) to evolve the model from z = 49 to z = 2. 



Instead of the u ^ UV background spectrum adopted by HKWM, we use the spectrum of 
Haardt & Madau (1996; hereafter HM] ), which is computed as a function of redshift based on the 



UV output of observed quasars and reprocessing by the observed Lya forest. The spectral shape is 
significantly different from v^^, but the UV background influences our Lya forest results primarily 



through the HI photoionization rate F, a cross-section weighted integral of the spectrum ( KWH , 
equation 29). The new simulation also includes star formation and feedback (see KWH), but this 
has no noticeable effect on the Lya forest results. A comparison between the galaxy populations 
of this simulation and the HKWM simulation appears in Weinberg, Hernquist, & Katz (1996). 

The mean opacity of the Lya forest depends on the parameter combination ^l/T. Since 
observational determinations of 0^ and T remain quite uncertain, we treat the overall intensity 
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of the UV background as a free parameter and scale it to match the observed mean Lya optical 
depth Tq,. When evolving the simulation, we divide HM's intensities by a factor of two, retaining 
their redshift history and spectral shape. We find that we must reduce the intensities by further 
factors of 1.28 and 1.38 at z = 2 and z = 3, respectively, in order to match the estimate 
fa = 0.0037(1 + of Press, Rybicki, & Schneider (1993; hereafter [PRSl). Although we apply 



this final reduction only at the analysis stage, to compute neutral fractions when generating 
spectra, the result is virtually identical to that of changing the intensity during dynamical 



evolution (Miralda-Escude et al. 1997, hereafter MWHK ) . In order to match the PRS mean optical 



depth at z = 3 with the original HM background intensity we would need rif, ~ 0.08, closer to 



the value advocated by Tytler, Fan, & Buries (1996; but see [Rugers Sz Hogan 199q and references 
therein). The value of fa plays the role of a normalizing constraint, used to fix the important 
combination of free parameters in our IGM model. Once ^i^/F is set, there is no further freedom 
to adjust the simulation predictions, and the remaining properties of the Lya forest provide tests 
of the cosmological scenario itself. 

We generate artificial spectra at z = 2 and z = 3 along 300 random lines of sight through 
the simulation cube, using the methods described in HKWM and Cen et al. (1994). We do not 
consider higher redshifts here because the strong absorption leads to severe blending of lines. The 
wavelength spread across the box is 23.4 A at z = 2 and 35.9 A at z=3. Each artificial spectrum 
contains 1000 pixels; an individual pixel has a velocity width ~ 2 kms^^ and a spatial extent 
20 comoving kpc, twice the gravitational softening length. In the Lya forest regime, the gas 
distribution is smooth on these scales. 



3. Fitting Voigt Profiles to Artificial Spectra 

We want the analysis of our simulated spectra to closely match that used in typical 
observational studies, HKCSR in particular. To this end, we have developed an automated 
Voigt-profile fitting routine, AUTOVP, which allows us to efficiently handle large quantities of 
simulated data and which provides an objective algorithm that can be applied to observational 
data. 

We add noise to our simulated spectra employing a combination of Gaussian photon noise 



with signal-to-noise ratio S/N = 50 in the continuum (corresponding roughly to the HKCSR data) 
and a fixed readout noise chosen to match the characteristics of the Keck HIRES spectrograph. 
Varying S/N changes our results only at the lowest column densities. While we know the true 
continuum level in the simulated spectra, this is not the case for the observational data. We 
therefore estimate the continuum in the simulated spectra by the iterative procedure commonly 
used for Echelle data: fitting a third-order polynomial to the data set, excluding any points lying 
^ 2(j below the fit, refitting the non-excluded points, and repeating until convergence is achieved. 
This technique is based upon the implicit assumption that the regions of lowest absorption in a 
high-resolution spectrum lie close to the true continuum level. Because the simulated spectra show 
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fluctuating Gunn-Peterson absorption that increases with z, continuum fitting has a systematic 
tendency to remove flux, an average of 1.2% at z = 2 and 5.7% at z = 3. The effect would probably 
be somewhat smaller in observational data because a typical HIRES Echelle order (~ 45A) is 
longer than one of our simulated spectra, giving a higher probability that the spectrum contains a 
region of genuinely low absorption. 

Given a normalized, continuum-fitted spectrum and its noise vector, we apply AUTOVP 
to detect lines and fit Voigt profiles. In its first phase, AUTOVP identifies lines and makes an 
initial estimate of their column densities and 6-parameters. Detection regions are identified above 



an 8a confidence level, following the method of [Lanzetta, Turnshek, fc Wolfe (1987)| . For line 



identification purposes, the data is convolved with a two-pixel-width Gaussian and la of noise is 
subtracted to yield a "minimum flux" . For non-saturated regions, a single Voigt profile is placed at 
the lowest flux value in the detection region, and A'^hi and b are reduced by small increments from 
large initial values until the model is everywhere above the minimum flux. For saturated regions a 
line is placed in the middle of the trough, and Nm and b are adjusted to fit the "cusp" regions, 
i.e. regions about five pixels wide on either side of the trough. The resulting first-guess line is then 
subtracted from the data to obtain a residual flux. Detections regions are identified in the residual 
flux, and the procedure is repeated until there are no more 8a detections. The line identification 
procedure is very robust, never failing for non-saturated regions and only occasionally producing 
a bad fit even in complex, blended, saturated regions, where A^hi and b are largely degenerate and 
the cusp regions are difficult to identify. 

In its second phase, AUTOVP takes the initial guess and performs a simultaneous x^- 
minimization on the parameters ("Ucentrai ) A^Hi , ^) of all lines within each detection region. Three 
independent minimization techniques are employed in conjunction in order to reliably identify the 
global minimum. AUTOVP then tries to remove any components with a formal error in A^hi 
or b comparable to the parameter value, refitting the detection region with one less line. If the 
resulting is lower than the original value the rejection is accepted, otherwise the fit returns to 
the original set of lines. AUTOVP thus attempts to fit the spectrum with as few lines as possible 
while still minimizing x^- If the fit after these line rejections is "good" (characterized empirically 
by X'^ ^ 2 per pixel), the program ends, otherwise it tries to add a line where the local contribution 
to x^ is greatest. The rare cases where AUTOVP fails to find a good fit are fiagged for possible 
manual intervention. AUTOVP is designed to interface with the PROFIT interactive Voigt-profile 
fitting package ( phurchill 1996| ). This graphical interactive fitter can be used to manually adjust 
poor fits, although this was required in so few cases that we base our simulation statistics entirely 
on the automated fits. 

In Figure || we show the results of AUTOVP applied to a S/N = 50, continuum-fitted 
spectrum at z = 3. The bottom panel shows the first-guess fit superimposed on the simulated 
spectrum. The top panel shows the final fit after the x^-minimization has been performed. 
Generally AUTOVP has greatest difficulty in obtaining first-guess fits in blended saturated regions 
like the one illustrated here. Nevertheless, the minimization produced an adequate fit with no 
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interactive adjustment. Lines with A'hi > 10 are indicated by the long tick marks above the 
spectrum, while lines with A'^hi < 10^'^ have short tick marks. The number of these small lines 
identified by AUTOVP is somewhat sensitive to the adopted S/N and the detection threshold. 
To keep our results fairly robust against details of our fitting procedure, we exclude these lines 
from our analysis and only focus on lines with A'hi > 10^^. AUTOVP has also been applied to 
observational data from H1216 and MgII2796, with results quite similar to those obtained from 
manual fitting. 

4. Results 

Figure |2| shows the column density distribution /(A^hi); the number of lines per unit redshift 
per linear interval of A^hi- Solid and dashed lines show the simulation results from AUTOVP 
at z = 3 and z = 2, respectively. The dotted line shows /(A^hi) obtained using the HKWM 
threshold algorithm at z = 3, with a fiux threshold of 0.7. As expected, the two methods yield 
similar results at high column densities, A^hi ^ 10^^'^ cm~^, but at lower A^hi AUTOVP deblends 
much more and finds many more lines. We find a similar trend at z = 2, though because of the 
reduced line crowding at lower redshift the agreement between the two methods extends down to 
A^Hi ~ 10^^ cm-2. 

Filled and open circles in Figure |2| show the observational results of PWRCL and 
HKCSR, respectively. The two determinations of /(A^hi) agree well in their regime of overlap, 
IQ13S ^^-2 ^ jy^^ ^ pj^-2_ rpj^^ j^j^^ g^jy resolution of the HIRES data allow 

HKCSR to detect much weaker absorption features, and their /(A'hi) continues to rise down to 
the lowest bin, A'hi ~ lO^'^'^ cm~^. Clearly a comparison between the simulations and HKCSR's 
published line statistics must be based on Voigt-profile fitting, since the vast majority of their 
lines lie in the region where line blending causes large differences between this method and the 
threshold algorithm. We compute the HKCSR /(A^hi) directly from their published line list, with 
no corrections for "incompleteness." HKCSR estimate such corrections from artificial spectra 
assuming an underlying model of randomly distributed, Voigt-profile lines with a power-law /(A^m) 
and a specified distribution of 6-parameters, and they conclude that their results are consistent 
with f{Ni{i) oc A^h/'^^ down to A^hi ^ 10^^'^ cm~^, where the correction for incompleteness is 
a factor of four. If we applied the same correction factors to the simulation results, the derived 
column density distributions would also rise in a nearly power-law fashion instead of turning over 
at low column densities. But the simulations provide no a priori reason to expect Voigt-profile 
lines or a power-law /{Nm), so we prefer to compare them directly to the data without trying to 
correct either for lines "lost" to blending. 

The mean redshift of the HKCSR lines is z = 2.9, so the closest comparison is to the z = 3 
simulation results. To make this comparison more exact, we convolved the z = 3 artificial 
spectra to a resolution AA = 0.06A (Av = 3.7 kms^^) before analysis, which has the minor 
eff'ect of removing some lines with A'^hi ^ 10^'^. When analyzed by Voigt-profile decomposition. 
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the simulation reproduces the large number of weak lines found in the HIRES spectra. In fact, 
the simulation overproduces the number of lines by a factor of 1.5 — 2 in the column density 
range 10^^ cm~^ ^ A^hi ^ 10^^ cm~^, a discrepancy that we will return to in §5. The rolloff 
at low column densities is also somewhat different, but the results for the weakest lines are the 
most sensitive to the details of the fitting procedure and to the modeling of noise and spectral 
resolution, so we regard this difference as less significant. 

In the regime where line blending is unimportant, /(A^hi) of the simulations drops by a factor 
of ^ 1.5 — 2 between z = 3 and z = 2. At low column densities /{Nyh) actually increases because 
of the reduced effects of line blending. As discussed in HKWM, Miralda-Escude et al. (1996), and 



MWHK , the evolution of the line population over this redshift range is driven primarily by the 
expansion of the universe, which lowers the physical gas densities in the absorbing systems and 
thereby lowers their neutral fractions and corresponding HI column densities. It is therefore more 
physically appropriate to think of /(A^hi) as evolving to the left rather than evolving downwards, 
though the quantitative effect is the same to within the accuracy of this simplified account. 

Figure |3| shows the distribution of 6-parameters for lines with A^hi ^ 10^^ cm~^ from HKCSR 
(solid histogram) and from the AUTOVP analyses of the simulation at z = 3 and z = 2 (solid 
and dashed curves, respectively). The 10^^ cm~^ cutoff eliminates lines whose identification and 
derived properties are sensitive to the value of 5"/ A^ or to details of the fitting procedure, though 
the results do not change qualitatively if we lower this cutoff to 10^^'^ cm~^. The threshold 
method (dotted curve) yields much larger 6-parameters than AUTOVP at z = 3 because many 
of its identified "lines" are extended absorption regions, which AUTOVP separates into narrower 
components. Table || lists the median, mean, and la dispersion of the 6-parameter histograms. 
The z = 3 simulation values for all three numbers are slightly larger than the HKCSR values, 
but the agreement is quite good given that the analysis procedures are not identical in all their 
details. The most significant difference in the distributions is the presence of many more narrow 
{b < 20 kms~^) lines in the simulation than in the data, a discrepancy that we discuss further 
below. 



5. Discussion 

Our most important result is that the CDM simulation reproduces the large number of weak 
lines found in HIRES spectra when it is analyzed by Voigt-profile decomposition. However, in the 
column density range 10^^ cm~^ ^ A'hi ^ 10^^ cm~^, the density of lines in the simulation at 
z = 3 is a factor of 1.5-2 higher than found by HKCSR at z = 2.9. Our simulation suffers from 
the inevitable limitation of finite numerical resolution, but the Lya absorbers are usually large, 
smooth, low-overdensity structures, and we would in any case expect higher numerical resolution 
to increase the number of lines rather than decrease it. This excess of lines may therefore indicate 
a failure of the 17 = 1, erg = 0.7 CDM model, a tendency to produce too much small scale clumping 
at z = 3. 
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An alternative possibility, quite plausible at present, is that we have set the intensity of 
the UV background too low given our adopted value of As discussed in §2, we choose the 
background intensity in order to match the PRS determination of the mean Lya optical depth, 
fa = 0.45 at z = 3. The statistical uncertainties in this determination are fairly small, but there 
are systematic uncertainties in the required extrapolation of the quasar continuum into the Lya 
forest region. We can match the HKCSR results if we increase the UV background intensity by 
a factor ~ 1.7, thus lowering HI column densities by a similar factor and shifting the simulation 
result for /(A'hi) to the left. This increase lowers the mean optical depth to fa = 0.32, which 
is well outside the la range of PRS (figure 4) at the HKCRS mean redshift z = 2.9 but is 
consistent with the PRS value at z = 2.65. It is somewhat above the value faiz = 3) ~ 0.25 
found by |Zuo &: Lu (1993) , who use a different data set and a different method of determining the 



quasar continuum. The uncertainty of our conclusions highlights the need for better observational 
determinations of faiz), which plays a crucial role in normalizing the predictions of cosmological 
simulations. The mean optical depth depends on the distribution of 6-parameters and on the small 
scale clustering of absorbers in addition to the column density distribution itself, so if fa is well 
known then the amplitude of /(A^hi) becomes an important independent test of the high-redshift 
structure predicted by a cosmological model. 

A second important result of our comparison is that the O = 1, us = 0.7 CDM model produces 
Lya forest lines with typical 6-parameters close to observed values. However, the simulation yields 
many more lines with 6 < 20 kms^^ than are found by HKCSR. A thermally broadened, single- 
temperature gas cloud produces a Lya absorption line of width b = 20(T/24, 000 K)^/^ kms~^. 
Low-6 lines arise in the simulated spectra because much of the absorbing gas is at temperatures 
of 10^ K or less, with its temperature determined by the balance between photoionization 
heating and adiabatic cooling ( 1VIWHK| ) . We will need tests at higher resolution to check that 



these temperatures are not artificially low because the simulation misses entropy production in 
unresolved shocks, but because such shocks would have to be quite weak, we expect that the effect 
of missing them would be small. 

Miralda-Escude & Rees (1994) pointed out that the process of reionization can heat the IGM 
significantly if it occurs fast enough to prevent radiative cooling losses. Our equilibrium treatment 
of photoionization (KWH) implicitly suppresses this effect, and since circumstantial evidence 



suggests that Hell reionization may have occurred at z ^ 3 (|Songaila &: Cowie 1996|) , we could 
be underestimating the gas temperatures at this redshift. The dot-dash line of Figure ^ shows 
the 6-parameter distribution obtained at z = 3 after adding 15,000 K to the temperatures of the 
SPH particles (a thermal energy equivalent to 4 Rydbergs per Hell photoelectron) , reducing the 
UV background intensity by 2.46 to restore fa = 0.45, then reextracting and reanalyzing spectra. 
Heating the gas eliminates the excess of low 6-parameter lines, though it worsens the agreement 
with HKCSR at 6 > 40 kms~^. We will investigate other treatments of reionization heating in 



future work, though the effects will be difficult to pin down because they depend on uncertain 
details of reionization ( Miralda-Escude fc Rees 1994| ). Model predictions for 6-parameters should 
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be more robust at z = 2, since by this time much of the energy absorbed during Hell reionization 
at z ^ 3 will have been lost to adiabatic cooling. 

Sharper tests of cosmological models against the statistics of the Lya forest can be obtained 
by expanding the redshift range of comparisons, by improving the determination of fa{z), and by 
applying AUTOVP to observational data, so that the analyses of simulated and observed spectra 
are identical in detail. We will also test models of the Lya forest using alternative statistical 
measures to characterize spectra, for if the physical scenario that emerges from cosmological 
simulations is correct, then Voigt-profile decomposition provides at best a rough guide to the 
density and temperature profiles of the absorbing gas. At the high S/N and resolution of the 
HIRES data, an absorption feature with a single flux minimum often shows asymmetries or broad 
wings, requiring two or more Voigt-profile lines to provide an adequate fit ( |HKCSR| ). Pairs of lines 
with small velocity separations have strongly anti-correlated 6-parameters, suggesting that many 
of these decompositions are not genuine physical blends ( Ranch 1996| ). While results such as 



these can always be accommodated within a discrete "cloud" model by postulating just the right 
clustering properties, they more likely signify the breakdown of the Voigt-profile paradigm itself, 
revealing the origin of the Lya forest in the diffuse, undulating gas distribution of the high-redshift 
universe. 
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Fig. 1. — Example AUTOVP fit (smooth curves) for a continuum-normalized, z = 3, artificial 

spectrum (histograms). The bottom panel shows the first-guess fit, while the top panel shows the 
final fit after x^-minimization. Long tick marks indicate lines with iVui > 10^^, while short tick 
marks indicate lines with iVni < 10^^. 

Fig. 2. — Column-density distributions /(A^m) = (PN/dz dNm, the number of lines per unit 
redshift per linear interval of HI column density. Solid and dashed lines show the AUTOVP results 
for 300 spectra along random lines of sight through the simulation at z = 3 and z = 2, respectively. 
The dotted line shows /(iVni) obtained by the threshold method at z = 3. Filled and open circles 
show the observational results of PWRCL and HKCSR, respectively. Error bars denote la Poisson 
counting errors. 

Fig. 3. — Distributions of 6-parameters, for lines with A?hi > lO^^cm"^. Solid and dashed curves 
show the AUTOVP results from the simulation at z = 3 and z = 2, respectively. The dotted line 
shows the result of the threshold method at z = 3; blending of features leads to very large line 
widths. The dot-dash line shows the AUTOVP result at z = 3 after gas temperatures have been 
increased by 15, 000 K. The solid histogram shows the HKCSR data. All distributions are computed 
as histograms with the same bins, but the simulation results are shown as curves to prevent visual 
confusion. 
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Tablc 1. Moments of ?)-paramctcr distributions. 



Spectra 


Algorithm 


z 


&med 


fcmean 


Ob 


HKCSR 




2.9 


31.4 


35.8 


16.3 


CDM 


AUTOVP 


3 


34.6 


39.3 


20.8 


CDM 


AUTOVP 


2 


30.6 


34.7 


16.5 


CDM+15K 


AUTOVP 


3 


40.9 


43.8 


19.5 


CDM 


Threshold 


3 


74.5 


83.9 


47.6 


CDM 


Threshold 


2 


43.4 


49.9 


26.3 
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